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Relativistic peculiarities at stepped surfaces: 
surprising energetics and unexpected diffusion patterns 
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We revive intriguing, yet still unexplained, experimental results of Ehrlich and coworkers 
[Phys. Rev. Lett. 77, 1334 (1996) and Phys. Rev. Lett. 67, 2509 (1991)] who have observed, 
that 5d adatoms distributed on (111) surface islands of 5d metals favor the adsorption at the clus- 
ter's edge rather than at the cluster's interior, which lies in contrast with the behavior of J^d and 
3d elements. Our state of the art ab initio calculations demonstrate that such behavior is a direct 
consequence of the relativity of 5d metals. 

PACS numbers: 68. 43. Be, 68.43.Jk, 68.65.-k, 73.20.At 



The recent surge of attention to nanoscale noble metal 
systems, serving as a base for nanostructured materi- 
als and paving the road to novel electronic devices and 
nanocatalytic agents, has reignitcd the community's in- 
terest in numerous peculiar properties of 5d elements. 
Among the abundance of experiments done over the past 
two decades on 5d metals some still remain not com- 
pletely understood despite fascinating and potentially 
fundamental results that have been obtained. 

In the present paper we revive the results of Ehrlich 
and coworkers^ who have observed, that the energetics 
and diffusion patterns of Pt on Pt(lll) and Ir or W 
on Ir(lll) display quite surprising features. In a clas- 
sical adsorption picture, known to be true for 3d and 4d 
metals 2,3 , single adatoms diffusing on stepped or island- 
covered surfaces at low temperatures tend to avoid areas 
close to descending step edges. Such repulsion is usu- 
ally attributed to the consequences of the electron den- 
sity smearing at step edges, known as the Smoluchowski 
effect. A sharp step in the atomic potential causes elec- 
trons to redistribute (flow down from the step) so as to 
smoothen out the electron density profile. The result- 
ing electronic depletion in the vicinity of a descending 
step edge makes adsorption there unfavorable, thus in- 
troducing a barrier (a few A wide) for diffusing adatoms. 
The depletion is strongest directly at the step so that dif- 
fusing atoms, once able to overcome the edge repulsion, 
tend to jump down onto the lower terrace or incorporate 
into the step^£. Such diffusion behavior can be described 
by a potential energy profile schematically presented in 
Fig. HJi. 

The first surprising feature found in experiments by 
Ehrlich and coworkers^ was that while the step repul- 
sion barrier was still observed for the diffusion of Pt on 
Pt(lll) and Ir or W on Ir(lll) its width was significantly 
increased (up to 10 A). Yet the most surprising find- 
ing, was the fact that once the temperature was raised 
high enough for the atoms to overcome the barrier intro- 
duced by the electronic depletion at the cluster's edge, 
the atoms did neither jump onto the lower terrace nor 
incorporate into the step but tended to migrate to the 
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FIG. 1. (color online) Schematics of the potential-energy 
landscape (a) which govern diffusion in 3d and 4d metallic 
systems^ and a landscape (b) which corresponds to the dif- 
fusion behavior observed in experiments on Pt and Ir 1 . 



step and diffuse along its upper edge without ever com- 
ing back to the cluster's interior. A potential energy 
profile which might result in such a diffusion pattern is 
sketched in Fig. [TJd (e.g. seeA). The novelty of the ob- 
served behavior has provoked a series of further studies 
on the subject^ and it has been suggested that the ob- 
served effects might be a consequence of unusual relax- 
ation patterns. Furthermore, similar phenomena have 
been observed in other surface systems. E.g., steps on Pt 
surfaces were found to facilitate the dissociation of 0?£ 
and CO 8 molecules and even promote water adsorption 
in step regions^. Yet a definitive answer as to the nature 
of such curious behavior has not been given. 

Here we demonstrate, that (i) the favorable edge ad- 
sorption can only be fully accounted for by taking into 
consideration the relativity of 5d elements (ii) but the 
width of the step-edge repulsion barrier on (111) surfaces 
of 5d metals is indeed dependent on atomic relaxations 
in the system. As a model system we consider Pt on 
Pt(lll) and generalize the obtained results to include 
miscellaneous other 5d and s-p species. 

To check whether atomic relaxations (as proposed ir>i) 
could indeed account for the behavior described above 
and at the same time have a flexible approach to rela- 
tivity we have made use of two different ab initio ap- 
proaches. To trace the effect of relaxations we have re- 
lied on the Vienna ab initio simulation package (VASP 
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where for geometry optimization a criterion of force-on- 
nuclei convergence to within 0.01 eV/A was used. We 
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FIG. 2. (color online) System setup for VASP (a) and KKR 
(b) calculations, (a) A stripe 14 atomic rows wide aligned 
along the [TlO] direction of a Pt(lll) surface, (b) Hexagonal 
Pt islands of ~ 25 A in diameter adsorbed homoepitaxially 
on a clean Pt(lll) surface, a single atom adsorbed along the 
shorter diameter of the island ([TT2] direction). 



have treated exchange and correlation effects with the 
local density approximation^ (LDA) and have used pro- 
jector augmented wave (PAW) potentials^. The details 
of the calculational routine can be found in Refs . 13 ' 14 . 
To model the experimental geometry the following setup 
has been used in our calculations. As the simulation 
of real-size islands in VASP is hardly possible from the 
computational point of view we have approximated ex- 
perimentally observed Pt islands^ by an infinite stripe 
14 atomic rows wide aligned along the [llO] direction of 
a Pt(lll) surface (Fig. [2k). Single adatoms have been 
adsorbed atop the stripe along the [112] line spanning 
the stripe from one step (type A) to the other (type B). 
Both fully relaxed and unrelaxed geometries have been 
considered. 

To obtain a better control over relativistic corrections 
we have employed an ab initio Korringa-Kohn-Rostoker 
(KKR) Green's function method in the atomic spheres 
approximation^. This method gives a proper description 
of physics on surfaces with atomic steps^t and islands^!. 
Normally the calculations are carried out under the inclu- 
sion of scalar-relativistic corrections to the Kohn-Sham 
equations. By excluding these corrections and carrying 
out fully self consistent calculations we were able to com- 
pare our system to its non- relativistic counterpart, thus 
underlining the governing role of relativity in the effects 
we describe. For KKR calculations we have constructed 
hexagonal Pt islands of ~ 25 A in diameter adsorbed ho- 
moepitaxially on a clean Pt(lll) surface (Fig. [2p). To 
study adsorption energies, single atoms have then been 
adsorbed along the shorter diameter of the island ([112] 
direction) . 

We start our investigation of the unusual energetics 
found in experiments*^ by attempting to reproduce ex- 
perimental observations in our VASP calculations. Let 
us first consider the adsorption of a Pt adatom diffus- 
ing across a Pt stripe (Fig. [2k). Corresponding ener- 
gies are presented in Fig. [3] for a fully relaxed (black 
squares) and ideal (red circles) geometries. As the en- 
ergy is determined up to an additive constant we present 
the energies respective to the adsorption energy in the 
center of the stripe in the relaxed configuration. The 
unrelaxed values have been further offset by 350 meV 
for clarity. Upon studying both relaxed and unrelaxed 
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FIG. 3. (color online) Adsorption energies of a Pt adatom 
diffusing in a straight line across a Pt(lll) stripe for a fully 
relaxed (black squares) and ideal (red circles) geometries. 
Atom position is given in numbered adsorption sites and 
A with respect to the center of the stripe. Integer numbers 
(filled squares and circles) denote fee and non-integer numbers 
(empty squares and circles) correspond to hep sites. 



configurations it becomes evident, that there is, indeed, 
an energetically unfavorable zone between the center of 
the cluster and its exterior. For the relaxed configura- 
tion the height of the barrier an adatom has to over- 
come to migrate from the center to the edge is of the 
order of 150 meV. Furthermore, in complete accord with 
the experiment 1,6 , adsorption at the edge of the island is 
more favorable by about 50 meV than that at the cen- 
ter. Comparing the energy profiles for ideal and relaxed 
geometries one notices that the presence of structural re- 
laxations strengthens the center-edge barrier effectively 
increasing its width. The results of KKR calculation fur- 
ther support the above-mentioned trends by yielding an 
adsorption energy at the B-step edge that is by 120 meV 
more favorable than the adsorption energy at the island's 
interior (the respective value for the A-step is 90 meV). 
If we consider the obtained energy differences in terms 
of Boltzmann statistics, adsorption at the step edge (at 
experimental temperatures of 100 K— ) is by a factor of 
300 more probable, than that at the center of a stripe or 
a cluster—. 

By studying the energy profiles in Fig. [3] we come to an 
evident conclusion that while atomic relaxations can play 
a role in defining the width of the edge repulsion barrier, 
which supports the hypothesis of Ehrlich and coworkers 1 , 
they can in no way account for the favorable adsorption 
at the step edge. 

Having ruled out the possibility of relaxations being 
the cause of preferable edge adsorption we have to ask 
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FIG. 4. (color online) LDOS at a single Pt adatom residing 
at the center of a hexagonal Pt(lll) island calculated non- 
relativistically (thin filled curves) and with scalar-relativistic 
corrections (thick solid curves) for s (red solid curves) and d 
(black broken curves) electrons. The s electron densities are 
scaled up by a factor of 20 for clarity. 

ourselves is what else distinguishes 5d metal steps from 
steps of other metallic systems. Indeed the remaining 
most obvious difference is the high mass of 5d elements, 
which makes relativistic effects most pronounced, and is 
known to cause interesting surface related effects 1 ^. 

The importance of relativistic effects for heavy-element 
chemistry and physics has been extensively discussed by 
Pyykko and Desclaux 2 - . In terms of electronic structure 
rclativi breaks down into several major effect a 20 ' 21 : (i) 
Electrons of s (and partially p) character experience the 
relativistic mass effect which causes the respective shells 
to contract^. For outer s and p shells the effect is some- 
what indirect as the relativistic contraction is only ex- 
perienced by the inner part of the wavefunction and the 
outer tails are pulled in as a consequence; (ii) palpable 
core levels' contraction leads to a stronger screening of 
the core, which decreases the pull on the d and / atomic 
orbitals causing them to expand radially and destabilize 
energetically 2 ^; (iii) the last effect is the spin-orbit cou- 
pling (SOC), which lifts the degeneracy of electrons with 
the same quantum number I. The SOC-related energy 
splitting may reach up to several me V 20 ' 22 , yet this is 
still at least an order of magnitude smaller than the en- 
ergy differences which we are trying to explain, so we 
will neglect the SOC for the rest of the present study. 
Effects (i) and (ii) can be accounted for by using scalar- 
relativistic corrections. 

Let us now see whether we can apply our knowledge 
about relativistic effects to explain the favorable step- 
edge adsorption on Pt islands. First of all we have to 
check what changes (in terms of the electronic structure) 
relativity introduces into our system. Let us consider 
a single Pt atom adsorbed at the center of a hexagonal 
Pt(lll) island some 25 A in diameter. To separate the ef- 
fects of relativity we perform KKR calculations both with 
scalar-relativistic corrections included and without them. 



The resulting local densities of s and d electron states 
(LDOS) for both cases are presented in Fig. 2) Thin 
filled curves represent the non-relativistic LDOS while 
thick solid curves reflect for the relativistic case (in both 
cases red solid curves denote the density of s electrons 
scaled up by a factor of 20 for clarity, and black broken 
curves represent the d-LDOS). Straight comparison im- 
mediately reveals both the s-level contraction (by about 
80 meV for the level localized around —0.75 eV) and the 
self-consistent d-level expansion (by about 300 meV for 
the level localized at —1.0 eV in the relativistic case). In 
accord with our expectations, the effect of relativity on 
the outer s-shell is weaker than the indirect effect of the 
core screening on the d-level. 

Such energy shifts of s and d states bring them closer 
together thus increasing their overlap and consequently 
the s-d hybridization. Much stronger, relativistically en- 
hanced, s-d bonding can well account for the difference in 
adsorption behavior between the elements of 5d and 3d- 
4d ranges. Indeed it is known, that the strength of the 
s-d hybridization can well decide over the preferences of 
bonding geometry of both relativistic and non-relativistic 
elements^ 3 -. 

To understand the impact of relativity and check our 
hypothesis about the predominant role of the s-d hy- 
bridization in the adsorption of relativistic impurities on 
relativistic substrates, let us consider, how the LDOS and 
the adsorption energy evolve as we move our Pt atom ad- 
sorbed on a hexagonal Pt(lll) island from the island's 
center to the edge (Fig. [SJ. If we neglect the scalar- 
relativistic corrections the s and d LDOS of a Pt adatom 
at the center of the Pt island (red solid and black broken 
thin filled curves in Fig. [5^, respectively) remains almost 
unchanged as we move the atom to a B-step edge of the 
island (thick solid curves in Fig. [5J red solid for s-LDOS 
and black broken for <i-LDOS). The only tangible effect 
is a slight destabilization of the d level and a slight re- 
duction of the s-d overlap, both of which are expectable 
due to the reduced coordinatio n 23 ' 24 . As a result, ad- 
sorption at the B-step for non-relativistic Pt is 80 meV 
less favorable than at the center. If we, however, include 
the relativistic corrections in our calculations the result 
is strikingly different. Already at the center of the is- 
land the s and d electronic levels of the Pt adsorbate 
are pushed closer together by the relativistic effects thus 
increasing the s-d overlap. But as we move the atom to- 
wards the step we find the both s and d levels are yet 
again shifted which lands them at the same energy and 
effectively locks them together. As a result, the s-d hy- 
bridization is maximized reducing the adsorption energy 
of the adatom and making edge adsorption favorable. 

Now only one questions remains to be resolved: why 
do the electronic levels shift upon moving the adatom 
to the edge in the relativistic case and remain almost 
unchanged when the relativity is shut down. This differ- 
ence can be best understood if one considers the bond 
formation at the center of the island and at its edge self- 
consistently. Reduced coordination at the edge reduces 
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FIG. 5. (color online) (a) LDOS at a single Pt adatom residing 
at the center (thin filled curves, red solid for s and black 
broken for d electrons) and at the B-step edge (thick curves, 
red solid for s and black broken for d electrons) of a hexagonal 
Pt(lll) island, calculated non-relativistically. (b) The same 
LDOS calculated with scalar-relativistic corrections. The s 
electron densities are scaled up by a factor of 20 for clarity. 

the number of bonds weakening the overall binding of the 
adsorbate to the substrate while making each of the re- 
maining bonds a bit stronger as one would expect in the 
context of Pauling's rules^. In the non-relativistic case 
the outer d orbitals are stable and cannot profit from 
the slight increase of single bonds' strengths and thus 
can not counteract the weakening of the overall bonding 
strength and the reduction of the s and p electron den- 
sity due to the Smoluchowski effect. In the relativistic 
case, however, the d orbitals are destabilized by the core 
screening and thus can more actively participate in the 
bonding (both s-d and d-d). In this way the adatom is 
able to profit from the bond enhancement by increasing 
the s-d overlap with the neighboring atom's shells, which 
in turn, self-consistently, brings the s and d levels closer 
together until they "lock" at the same energy, producing 
the LDOS displayed in Fig. [5p. 

The feasibility of such an adsorption and bonding sce- 



nario is further supported by the fact, that mixed sys- 
tems, e.g. a 3d adatom on a 5d substrate, show an in- 
termediate behavior, i.e. the s and d orbitals are more 
susceptible to the adsorption geometry changes than in 
non-relativistic systems, yet the introduced shifts of s 
and d levels are smaller than e.g. in Pt/Pt(lll). For 
example, for Cu/Pt(lll) the s-d interlevel energy dif- 
ference is reduced by 25% upon the shift of the adatom 
from the center of the island to its edge. The s-d overlap 
is thus slightly increased. However, this increase turns 
out to be insufficient to afford edge adsorption: accord- 
ing to our calculations edge adsorption for Cu/Pt(lll) is 
by about 100 meV less favorable than the adsorption on 
a clean surface or in the center of an island. 

If we stray from our model system of Pt on Pt(lll) 
we would also notice, that similar behavior can be found 
in many systems where relativity is present or the outer 
orbitals of adsorbates are unstable enough to be strongly 
affected by the bonding strength increase at a step edge. 
For example, our VASP calculations show that while 3d 
and 4d elements (Cu, Au and Pd) tend to adsorb more 
readily in the interior of a Pt(lll) island, relativistic 5d 
species (Pt, Ir) and species of a more pronounced s-p 
character (C, O, N) follow suit of Pt and favor step edge 
adsorption. 

We might thus conclude that the curious diffusion pat- 
terns observed in experiments in Refs. -' 6 is a consequence 
of a delicate balance where on one side of the scales the 
Smoluchowski effect and coordination reduction strive to 
weaken the bonding of adsorbates to the substrate and on 
the other side the individual bonds' enhancement striv- 
ing to outweigh them. Now we have shown that while 
in non-relativistic system the edge adsorption is unfavor- 
able, relativistic effects are well able to topple the bal- 
ance. 

Apart of the fundamental value of this notion, we be- 
lieve that it might prove to be applicable in surface engi- 
neering, e.g. for the bottom- up creation of various surface 
based nanostructures for further electronic and nanocat- 
alytic research^. 
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